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ABSTRACT 


A  complete  analytical  description  of  a  device  to  measure  angle 
of  arrival  is  presented  along  with  predictions  of  expected  results 
for  various  input  apertures  and  turbulence  conditions.  This  work  is 
necessary  for  defin'ng  experiments  made  with  the  instrument  of  RADC. 
Also  presented  are  some  extensions  of  previous  work  on  temporal 
spectra  of  phase  difference  measurements.  Specifically  the  case  for 
slanted  paths  is  considered. 
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INTRODUCTION 


This  is  tne  second  quarterly  technical  report  under  Contract 
F30602-72-C-0305.  The  object  of  the  contract  is  to  provide  theoretical 
support  to  the  experimental  program  in  atmospheric  propagation  of 
10.6  micron  light  beams  at  the  Rome  Air  Development  Center.  Specific 
problems  such  as  theoretical  prediction  of  temporal  spectra  of  light 
beam  amplitude  and  phase  fluctuations,  averaging  time  required  for 
optical  atmospheric  data  have  been  the  subjects  of  past  investi¬ 
gations.  The  work  to  be  reported  on  this  quarter  is  a  theoretical 
investigation  of  the  apparatus  used  for  the  measurement  of  angle  of 
arrival  at  RADC.  A  complete  analytical  description  of  device  operation 
is  presented  along  with  predictions  of  expected  results  for  various 
input  apertures  and  turbulence  conditions.  This  work  is  necessary  for 
defining  experiments  and  interpreting  results  of  measurements  made  with 
the  instrument. 

Also  presented  are  some  e  .ensions  of  previous  work  on  temporal 
spectra  of  phase  difference  measurements.  Specifically  the  case  for 
slanted  paths  is  considered.  This  work  applies  to  measurements  made 
when  the  optical  receiver  is  in  the  receiving  tower  at  the  RADC  PATS 
site,  and  for  vertical  beams. 

The  angle  of  arrival  apparatus  will  be  described  first. 


ANGLE  OF  ARRIVAL  APPARATUS 

The  following  is  a  description  of  the  apparatus  used  at  RADC  to 
measure  angle  of  arrival.  The  description  was  started  in  conjunction 
with  a  report  on  angle  of  arrival  (3163-1)  put  out  under  a  previous 
contract.  The  description  was  put  aside  to  be  completed  at  a  later 
date,  a  task  which  is  performed  herein. 

The  description  was  undertaken  to  provide  a  better  understanding 
of  an  arrival  angle  angle  meaning  device  designed  and  constructed  for 
RADC  by  the  Perkin  Elmer  Company.  Descriptions  of  the  device  given  in 
their  literature  (PEI)  are  generally  qualitative.  In  order  for  the 
measurements  to  be  meaningfully  compared  with  theoretical  predictions 
appearing  in  the  literature  (Hufnagel ,  1963),  (Fried,  1966),  (Heidbredder, 
1967),  (Tatarski ,  1962)  and  to  obtain  a  more  quantitative  understanding 
of  the  device  operation  this  analysis  was  undertaken.  Indeed  as  a  result 
of  this  work  we  find  that  the  analyticil  expression  for  the  quantity 
measured  is  not  the  same  as  any  of  those  commonly  chosen  to  define 
angle  of  arrival.  Furthermore  the  weighting  of  light  entering  the 
input  aperture  is  also  found  to  be  quite  unique. 
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Angle  of  arrival  measurements  are  in  general  useful  because  they 
give  angular  coordinates  of  objects.  One  generally  defines  angle  of 
arrival  as  the  angles  defining  the  normal  to  the  wavefront  of  a  light 
wave  originating  from  a  point  object.  If  there  are  no  intervening 
index  fluctuations  then  the  wavefront  will  be  a  spherical  one  centered 
on  the  object  with  normals  indicating  the  direction  to  the  < bject.  If 
there  are  interviening  refractive  index  fluctuations,  then  the  arriving 
wavefront  may  be  different  from  spherical  or  may  be  spherical  with  a 
different  center  or  radius  of  curvature,  giving  a  false  measurement. 

The  device  senses  angle  of  arrival  fluctuations  by  examining  the 
position  of  the  focused  diffraction  limited  point  image.  The  trans¬ 
verse  displacement,  s,  is  measured  and  results  are  interpreted  in  terms 
of  angle  of  arrival,  a,  by  the  standard  relationship  a  =  6/c  where  d 
is  the  telescope  image  distance. 

All  angle  of  arrival  measurements  are  limited.  Even  in  the 
ideal  case  measurements  having  better  accuracy  than  the  classical 
angular  resolution  determined  by  the  input  aperture  of  the  optical 
system,  1.22  x/D  for  a  round  aperture,  are  difficult.  If  there  are  no 
refractive  index  fluctuations  the  object  of  a  stationary  point  source 
is  a  stationary  diffraction  pattern  the  position  of  the  object  being 
determined  by  the  center  of  the  pattern.  If  there  is  a  spherical  wave- 
front  with  center  modified  by  the  index  fluctuations,  then  the  diffraction 
pattern  will  be  cenered  at  a  different  spot,  the  amount  of  displacement 
being  proportional  to  the  apparent  angular  deviation  of  the  object.  If 
there  is  a  nonspherical  wavefront,  then  the  standard  diffraction  limited 
pattern  is  broken  up,  the  amount  of  degradation  depending  on  the 
aperture  diameter  compared  with  the  mean  transverse  size  of  the  atmospheric 
fluctuations.  The  device  determines  the  mean  spot  position  and  averages 
over  the  spot  shape. 

The  apparatus  is  shown  schematically  in  Fig.  1.  Light  from  a 
point  source  degraded  by  the  turbulent  atmosphere  is  collected  by 
the  receiving  telescope  and  Imaged  to  a  degraded  diffraction  limited 
spot.  The  spot  is  then  chopped  with  the  scanning  reticle  wheel  and 
applied  to  a  detector.  The  detector  output  has  the  form  of  a  periodic 
signal  as  shown  in  Fig.  2.  The  scanning  rate  of  the  wheel  is 
sufficiently  fast  so  that  the  spot  is  chopped  many  times  before  the 
atmosphere  changes  it  appreciably.  The  detector  output  is  then 
filtered  to  pass  the  fundamental  frequency.  This  filtered  signal  is 
compared  for  phase  with  a  standard  signal  derived  from  the  opposite  side 
of  the  same  scanning  wheel.  The  phase  difference  is  normalized  to  zero 
mean,  squared  and  averaged.  The  mean  square  phase  is  then  related  to  the 
angle  of  arrival  variance  in  the  following  analysis. 

At  some  instant  of  time  the  voltage  output,  S,  of  the  signal 
detector  will  be  given  by 
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Fig.  1.  Angle  of  arrival  apparatus. 
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dx  |Eo(p)| 


p  denotes  a  point  (x,y)  in  the  image  plane, 

EQ(p)  =  light  field  in  the  image  plane 

CQ  =  constant  relating  field  to  output  voltage.  C0 

contains  such  factors  as  detector  quantum  efficiency, 
impedance  of  free  space,  output  impedance,  etc. 


The  slits  on  the  reticle  wheel  extend  sufficiently  far  in  the 
ydirection  so  as  to  contain  the  entire  spot.  The  integration  limits 
x  and  x  denote  the  x  values  of  the  edge  of  the  slit  in  front  of  the 
spot  at  the  instant  under  consideration. 


Using  the  Fourier  transform  relationship  between  the  field, 
El(r),  in  the  telescope  aperture  plane  and  the  image  plane  field’  the 
expression  for  the  signal  becomes 


<2)  s  =  co(^)  LdyLdxlL 


,  00  _  _  _  _  _  jip  •  (r,  -r0 ) 

drl|  dr2W(rl  )W(r2)El  (rl  )Ef^r2)e  d  2 


where 


r-|  and  r2  denote  positions  (rlx,rly)  and  (r2x,r?v)  in  the 
aperture  plane.  y 

W(r)  is  the  aperture  shape  function  W=0  outside  the 
aperture  and  W=1  inside  it. 

A  =  2ir/k  is  the  light  wavelength 

d  =  receiving  lens  image  distance. 

The  first  step  is  the  performance  of  the  x  and  y  integrations. 
These  are  straight  forward,  the  results  being 
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•  H(r,)W(r2)E(r,)E»(r2). 
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Equation  (3)  is  put  in  more  tractable  form  b>  „he  transformation  to 
sum  and  difference  coordinates.,  both  in  the  pupil  plane  and  in  the  image 
plane.  Thus  we  define 


(4) 


-  x“-x' 


xH+x' 
o  -  -  — 


5  =  rrr2 

-  Vr2 

n 5 — 


6  is  the  scanning  slit  width  and  o  the  slit  center.  At  this 
point  it  is  also  appropriate  to  separate  out  a  flat  portion  from  the 
wavefront.  We  define  E-j  (r)  as  the  incoming  field  with  the  phase 
measured  from  a  plane  normal  to  the  vector  k. 


(5) 


E-j  (r)  =  E(r)  e 


r\ 


The  particular  choice  for  will  be  indicated  later.  (The  apparatus 
averages  over  <y.)  In  terms  of  these  variables,  Eq.  (3)  becomes 


(6) 


k  C  6 

s<‘>  ■ 
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-  ,  ,  ^d°+,cx^x+^KvCv  k^x6 

dn  0(r  )  e  J  x  x  y  ^inc^) 


W(n  +  §■)  W(n  -  f)  E^n  +  f-)  E*(n  -  f). 


6 (Cy )  function  indicates  that  the  apparatus  is  insensitive  to 
wavefront  tip  in  the  y  direction  and  the  spot  can  have  any  vertical 
position  in  the  slit. 

We  now  incorporate  the  slit  motion  into  the  analysis  by  letting 
the  slit  mean  position,  o,  move  linearly  in  time  with  constant  velocity, 
v,  across  the  spot. 

(7)  o  =  vt  . 

The  slit  is  assumed  to  be  wider  than  the  spot.  We  imagine  that  a 
humped  shaped  signal  is  produced  as  the  slit  moves  across  the  spot.  This 
gives  one  cycle  of  a  periodic  pattern  that  is  repeated  a  large  number 
cf  times.  The  wheel  velocity  is  assumed  sufficiently  fast  so  that  a 
lary  number  of  slits  scan  the  spot  before  the  spot  shape  changes 
appreciably,  thus  generating  an  essentially  periodic  function  of  period 

(8)  t  *  26/v. 
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The  periodic  signal  will  then  be  represented  in  the  form  of  a  Fourier 
series.  We  further  adjust  the  time  scale  so  that  the  signal  peak  will 
come  very  close  to  zero  time.  This  is  done  by  adding  in  a  time  bias 
so  that  the  coefficient  of  £x  in  Eq.  (6)  is  zero  for  zero  time.  To  do 
this  we  put 


(9) 


k  k  .  k 

jfl  •  K  =  -jVt  -  K  =  jVt 

d  x  d  x  d 


This  is  equivalent  to  delaying  the  time  scale  so  that  it  coincides  with 
a  spot  that  has  been  moved  sideways  by  an  amount  kx  d/k  due  to  the 
mean  plane  wavefront  tilt  coivonent,  kx. 


The  form  of  the  Fourier  series  to  be  used  is  that  in  Eq.  (10) 

k  d 

-  j*n  iijyt  p  j*  j^(f  +  £-) 

(10)  S(t)  =  l  an  e  n  e  e  =  J  an  e  n  e  6  kv 


where 


(ID  an  e 


__  eJ7rn 

J  -t>/\ 


k  d 

X 


ntrvt1 


j 

S(t')  e  0  dt 1 . 


The  limits  are  the  natural  ones,  i.e.,  those  centered  at  t'=0. 

Inserting  Eq.  (7)  into  Eq.  (6),  Eq.  (6)  into  Eq.  ('ll)  eliminating 
v  and  performing  the  time  integration  gives 


(12) 


an  e  2irH  e 


11  «  11 


—  .  J<v£v 

dn  5(Cy)  e  y  y 


•  W(n  +  |-)W(n  -  |-)E^(n  +  |-)E^(ri  -  |-)  s1nc(-2^-)sinc(^f v 

We  next  consider  the  function  performed  by  the  electronic  pro¬ 
cessing  equipment.  The  bandpass  filter  selects  out  the  signal 
component  at  the  fundamental  frequency.  The  phase  is  measuied  by  means 
of  zero  point  crossings.  The  phase  difference  from  the  reference 
signal  Is  found  and2reduced  to  zero  mean,  and  the  mean  square  phase  dif¬ 
ference,  <Un"4>ref/  >>  Is  obtained.  We  want  to  derive  an  expression  for 
that  mean  square  phase  difference  relating  It  In  a  simple  way  to  the  mean 
square  angle  of  arrival.  The  procedure  w*  1 1  be  to  develope  an  equivalent 
situation  and  show  how  that  leads  to  the  desired  results. 


-nir 


6 


The  first  step  is  to  note  that  the  measurement  of  the  zero 
crossings  would  be  unaffected  if  a  feedback  circuit  and  a  variable 
qain  amplifier  were  added  to  the  system  to  maintain  a  constant 
amplitude  for  the  fundamental  harmonic.  That  would  be  equivalent  to 
letting  the  constant  C0  be  a  random  function  defined  as  being  pro¬ 
portional  to  the  reciprocal  of  the  amplitude. 


We  also  choose  at  this  time  a  particular  definition  of  the  wave- 
front  normal  component,  kx.  Specifically,  <x  is  chosen  to  make  the 
imaginary  part  of  the  integral  in  Eq.  (12)  zero.  That  is,  denoting 
the  integral  in  Eq.  (12)  by  I ,  we  have 


03) 


Im  I  =  Inn 


d£ 

-  T  k6^X  .  i 

•  E-| (n  +  §-)  E*(n  -  §-)  sinc(-^)  sinc(-j - *)>■  -  0. 


dn  <5Uy)  e  y  y  W(n  +  f-)  W(n  -  |) 


Taking  the  magnitude  of  both  sides  of  Eq.  (11)  with  the  imaginary  part 
of  the  integral  equal  to  zero  then  gives 


(14) 


k  C  6 

-fr— — Re  ( I )  =  constant. 
2nd 


Imagine  next  that  we  detect  the  first  harmonic  signal  giving  as 

output  ai  eJ<f>1  and  take  a  long  time  average  of  it.  Assume  that  the 
average  is  sufficiently  long  so  as  to  have  the  ergodic  hypothesis 
hold,  and  we  equate  ensemble  and  time  averages.  Since  the  magnitudes 
are  constant,  the  result  would  be 


(15) 


ve 


jni 


<e 


Kxd 


k  C  6 
0 

2nd 


Re(I). 


Canceling  out  the  constant  factors,  assuming  a  zero  mean  and  normal 
distribution  for  the  arrival  angle  fluctuations,  and  using  the  re¬ 
lationship  <eju>  *  el/^<u  >  for  a  normally  distributed  zero  mean  random 
variable,  u,  gives 


,.c,  ,2  _  nd  2 

O6)  «t>i>  -  tt<Kx  * 

Equation  (16)  is  the  desired  expression  relating  the  phase  variance 
to  the  angle  of  arrival  variance.  Equation  (16)  is  exactly  what  one 
would  expect  on  simple  grounds.  A  given  value  of  <x  would  produce  a 
spot  sideways  shift  of  distance  d<x/k.  This  would  then  correspond  to 
an  electronic  phase  shift  of  <j>-j  =  (tt/ 6 ) (dicy/k ) ,  in  accord  with  Eq.  (16) 
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We  now  want  to  continue  with  an  examination  of  the  main  require¬ 
ment  in  obtaining  Eq.  (16),  that  the  imaginary  part  of  the  integral  in 
bq.  (12)  be  zero,  as  indicated  in  Eq.  (13).  We  replace  Fi(r)  in  terms 
of  its  original  expression  and  substitute  for  the  incoming  field  in  terms 
of  a  phase  and  log-amplitude 

(17)  E(F)  =  et(F)  +  j*(F). 

We  further  assume  that  the  major  part  of  the  phase  fluctuations  are 
contained  in  the  wavefront  tilt  (Fried,  1966)  so  that  the  difference 
between  the  flat  plane  and  the  actual  wavefront  is  less  than  a  tenth 
radian.  The  sine  of  the  phase  difference  (4>-KxX-Kyy)  can  then  be 
replaced  by  its  argument.  Equation  (13)  then  simplifies  to 


(18) 

where 

(19) 


'*  =  h/l 


N  = 


7-7  k65x  •  k65x 
dn«U  )W(n  +  |)W(n  -  ^sinct-^sincf-g - 


fc(n  +  §-)  +  *-(n  ■  t)  _  r  7X  x 

!  (<f>(n  +  y)  -  <t*(n  -  §-) ) 


(20) 


A  = 


—  _  —  kfi£  k6E 

dn«Uy)W(n  +  f)W(n  -  f-Jsinct-^sinct-g - it) 


Jt(n  +  f)  +  a(i  “  f") 


The  expressions  in  Eqs.  (19)  and  (20)  can  be  simplified.  N  is 
the  sum  of  two  integrals  corresponding  to  the  two  <p  functions  in_the 
last  parenthesis.  We  first  express_Eq.  (19J  in  terms  of  rj  and  r2 
using  Eq.  (4)  and  then  interchange  rj  and  rg  in  the  second  integral. 
All  the  factors  remain  identical  with  those  in  the  first  integral 
except  for 

^(r„i-rv9;  \ 

•*] 


sine 


P 


xl 

T 


which  goes  into 


/k6(rrr^)  \ 

sinc[ - +  *) 
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The  result  is  Eq.  (21).  The  expression  for  a  in  Eq.  (20)  can  be 
similarly  simplified  since  5x(=rxrrx2)  behaves  the  saire  as 


4> (n  +  §•)  -  4> ( n  ■  p)  (=4> ( ^] )  -  ♦  (^J) 


in  Eq.  (19)  under  exchange  of  r-|  and  r2.  The  result  is  Eq.  (22) 

t(r,)  _ 

.  .  I—  V  I  I 


(21)  N  - 


(22)  A 


where 


f  _  _  v » i  /  _  _ 

dr.|W(r.|)e  4>(r-j  )F-|  (r-j ) 


dr]W(r1)e  rxlFl^V 


(23)  >:-|  ( r-j )  =  }{dr2W(r2)6(ryl-ry2)e 

.  fk6(rxTrx2M 

x  sinc^ - 23 - H 


n(r2) 


i .  Jk5<rxrrx2> 

'sine - 3 - 


1  .  (kS‘rx1-rx2>  .  „ 

ir  -Sincl - ^ -  +  TTlr 


The  derivation  is  now  restricted  to  cases  where  amplitude 
fluctuations  of  the  incoming  beam  are  neglegible.  This  situation  is 
close  to  that  met  for  10.6  micron  light  at  ranges  found  at  the  RADC 
propagation  range.  With  this  restriction  the  coefficient,  a,  of  kx 
becomes  constant  over  the  ensemble  enabling  one  to  write  an  expression 
for  its  mean  square  value. 


(24) 


<v. 


>  = 


<n2> 

.2 


x.2 


a 2  as 


Our  final  job  is  now  to  evaluate  expressions  for  <N  >  and 
functions  of  geometry  r.nd  coherence  parameters.  Squaring  and  averaging 
the  expression  for  N  in  Eq.  (21),  substituting  for  the  correlation 
function  in  terms  of  the  relationship 


(25)  <♦(?,)♦(?,•)>■  yvv>  *  W*1  •  tyrrri,)) 

and  recognizing  that  the  integral  with  D^®)  has  value  zero  gives 
Eq.  (26) 
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(26) 


<N2> 


rdr. 


dr„ 


dri 


dr2<5  (  ryl  "ry2  }  5  (  ry  l"ry2  ) 


W(r-j  y W(r2) W(r-j  )W(r2)si nc(^g-(rxl “r^Jsi nc(^-( ~r*2) 


,k6 


(s1nc(^-(rx]-rx2)-,)-s1nc(^(rxrrx2)+,)i  x 


,k6 


'ts,'nc(r<r;)-r;2)-»)-s1nc(^(r;)-p^)+,)}  x  D^ryrj). 

Equation  (26)  can  now  be  partitioned  because  Mri-H )  is  independent 
of  r£  and  r£.  Separating  the  r2  and  r£  integrations  gives 

(26a)  <N2>  =  -  y  ~ 


dri  Jdri  Vrrri  )w(ri  )p(ri  Wrj ) 


where 


(26b)  F(r-| )  =  jdr26(ryl-ry2)W(r2)sinc(||(r 


r  k  6 


xrrx2,) 

k  6, 


IS 


•{sinc(5-(rxi-rx2)-’T)-sinc(d-(rxi-rx2)+Tr). 

p 

An  expression  for  A  can  be  similarly  derived  from  Eq.  (22).  It 


(26c)  a2  =  jdr1Jdr]rlxr]xW(r1)W(r')F(r1)F(F'). 

Equations  (25)  and  (26)  will  thus  be  evaluated  to  give  an  analytical 
expression  for  <k2>.  *  3 

^  .  Before  proceeding  with  the  numerical  expressions  for  <N2>  and  a 2 
it  might  be  interesting  to  compare  the  expression  thus  obtained  for 

comPrtr®x1e  expressions  appearing  in  the  literature  (Fried,  1966) 
(Heidbreder,  1967).  The  commonly  suggested  expression  for  «2>  is 

(27)  .  -  ?ffdfldFiD 

X  _ 

Jdr1dr|x2x]2W(r1)W(rj) 
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The  form  of  tne  two  expressions  is  indeed  identical,  the  difference 
being  the  replacement  of  the  odd  function  xi  by  another  odd  function 
F(r-| )  will  be  examined  further  later. 

We  now  proceed  with  numerical  evaluation  of  <N2>  and  a2  The 
first  step  is  the  nondimensionalization  of  the  relaverit  expressions 
with  coordinate  substitutions,  a  choice  of  (p ) , 

(28a)  N  =  k  6D/2d 

(28b)  u  =  (rxl-rx2)  N/2D 

(28c)  X1  =  2rxl/D 

(28d)  Y1  =  2ryl/D 

(28e)  R1  =  (X1,Y1) 

(28f )  D^p)  =  6.88(p/r0)5/3, 


and  determination  of  the  integration  limits.  The  results  are 


(29a)  <N2>  =  -  6.88(—  )5/3  £ 


r 


1  +J1-X 


dX| 

-1 


dY, 


dXJ 


'1 


RJ-i  -J i -xi2 


f  \  (Rt^  )5/3F(R1  )F(R1- ) 


(29b) 


F(R1) 


D 

N 


r^x,+JTvf) 

du  sine  u  (sinc(2u-Tr)  -  sinc(2u+ir)} 


(29c)  ^  -  +  D 


dX, 


-1 


•i  1-Xt  1 

'  r 

dY]  dX 


-ji- 


i 


+  l-Y 


72 

l 


dY1,X1X^F(FTJ)F(Rp. 

i' 

1 


- ? 

-  l-Y 
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It  is  interesting  to  note  that  the  form  resulting  from  the  nondimensional- 
ization  is  precisely  that  predicted  on  the  basis  of  very  simple  arguments 
multiplied  by  a  correction  factor  depending  on  the  size  of  the  aperture: 


(30) 


Kx  3.44 

~~2 —  =  TT 
k  k  u 


D»(D) 

2kV~ 


1 


The  nondimensional  comb  nation  N=k6D/2d  =  ttD/ ( x/6 )d  is  a  type  of  Fresnel 
number  composer1  of  (x/6)d,  that  is  the  diffraction  limited  length  in  the 
aperture  plane  due  to  the  slit  width  6  in  the  image  plane.  Then  D/ ( x/6 )d 
is  the  number  of  diffraction  1 i mi ti e  lengths  in  the  aperture  diameter 
or  roughly  the  number  of  independent  regions  across  the  diameter. 

An  interesting  point  about  the  apparatus  can  be  indicated  by 
examining  the  function  F(X,Y).  To  do  this  we  first  sketch  the  integrand. 
This  is  done  in  Fig.  3.  The  limits  on  the  integral  in  F(X,Y)  are 
symmetric  about  NX/2  and  extend  a  distance  (N/2){l-Y?  in  either  direction 
about  NX.  Examples  of  various  cases  for  the  limits  are  shown  in  Fig.  3. 
In  case  I  in  Fig.  3,  NX/2  is  fairly  small  and (N/2^1-Y^  is  large  so  that 
the  integration  limits  include  a  fair  portion  of  the  integrand. 

The  majority  of  the  area  is  under  the  two  central  lobes  so  that 
the  integral  with  limits  shown  in  case  I  will  be  quite  small,  because 
the  positive  and  negative  lobes  cancel.  However,  if  the  limits  are  as 
shown  in  cases  II  and  III  then  the  integral  will  be  larger  because  the 
region  of  integration  includes  only  one  of  the  central  lobes.  This  case 
occurs  when  the  semi  extent  of  the  integration  region,  N/2 Jl -Yz  nearly 
equals  N/2  X, 

(31)  Sr-=  jkl-Y2)172, 

or,  using  Eqs.  (28c)  and  (28d)  when 


That  is  F(X,Y)  is  large  when  the  integrand  in  the  integrals  for  N  and 
A  are  near  the  aperture  edge.  Indeed  using  the  fact  that  the  first 
2ero  in  the  integrand  of  F(X,Y)  occurs  for  u=tt  it  is  simple  to  show 
that  F  is  large  in  a  ring  within  a  distance  roughly  d(x/6)  around  the 
aperture  edge.  The  net  result  is  that  light  entering  the  center  of 
the  aperture  has  no  effect  on  the  average.  Indeed  the  apparatus 
might  operate  more  effectively  with  a  better  signal-to-noise  ratio  if 
an  obscure  form  were  inserted  to  block  the  aperture  center. 
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Fig.  3.  The  integration  integrand  and  three 
examples  of  the  integration  limits. 

The  integrals  in  <N^>  and  a2  can  further  be  evaluated  for  very 
small  and  very  large  N.  The  significant  features  lie  in  the  limiting 
expressions  for  F-j  (IT) .  For  example  for  small  N 


This  occurs  iri  both  <N2>  and  so  that  the  ratio  <N2>/a  is  independent 
of  N.  It  can  also  be  shown  that  for  large  N  <N2>/a2  is  similarly 
independent  of  N. 
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2  2 

The  integrals  in  <N  >  and  a  have  been  calculated  numerically 
for  the  case  of  large  aperture.  For  X  =  10y,  D  =  40  cm  and  6  =  N  = 
46.103,  f (N)  =  .975.  For  N  =  92.203  f(N)  =  .971  illustrating  the 
constancy  of  f (N )  for  large  N. 

To  summarize  this  section,  a  theoretical  description  of  the 
particular  angle  of  arrival  measuring  apparatus  has  been  presented 
for  the  case  of  small  amplitude  fluctuations.  Several  conclusions 
are  indicated. 


1)  The  device  does  not  measure  what  has  been  commonly  defined 
as  angle  of  arrival,  but  does  measure  something  with  similar 
aperture  symmetry  properties. 


2)  The  basic  dependence  of  arrival  angle  variance  on  aperture  D, 
wavelength,  27r/k,  and  coherence  length,  r0,  is  that  pre¬ 
dicted  on  the  basis  of  simple  arguments.  It  is 


<„2>  = 


0d(D) 

TT' 

\c\r 


3)  For  very  large  and  very  small  apertures,  arrival  angle 
variance  is  independent  of  N  =  ksD/2d  (6  =  scanning  wheel 
slot  size,  d  =  image  distance). 


4)  Only  light  entering  the  aperture  in  a  ring  of  width  dx/6  at 
the  aperture  edge  is  effective  in  determining  the  arrival 
angle  variance. 


PHASE  DIFFERENCE  TEMPORAL  SPECTRA 
ALONG  SLANT  PATHS 

The  temporal  spectrum  calculations  of  phase  difference  fluctuations 
for  the  case  when  the  wind  is  parallel  to  the  phase  point  separation  have 
been  extended  to  include  the  cases  of  vertical  and  slant  (upward  and 
downward)  paths.  Input  assumptions  include  turbulence  strength  C^  is 
proportional  to  altitude  to  the  -4/3  power  and  the  outer  scale  proportional 
to  the  altitude.  The  wind  is  assumed  to  be  constant  along  the  path. 

This  extends  work  of  previous  reports  (3432-1). 

The  calculations  have  been  separated  so  that  the  slant  path 
spectrum  is  obtained  from  the  horizontal  path  spectrum  for  a  path  of 
the  same  length  and  located  at  the  lower  end  point  of  the  slant  path. 

The  additional  factor  (h)  is 

m  h .  H."ila 
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where  Hgax  is  the  upper  end  pcrnt  altitude  and  Hm-jn  is  the  lower  end 
point  altitude. 

The  high,  intermediate  and  low  frequency  regions  of  the  horizontal 
path  spectrum  are  given  in  Table  I.  Note  that  for  large  normalized 
separation  there  are  only  the  high  and  low  frequency  regions. 


TABLE  I 


PHASE 


DIFFERENCE  SPECTRUM  WHEN  THE  WIND  DIRECTION  IS  PARALLEL 
TO  THE  LINE  CONNECTING  THE  PHASE  POINTS 
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Figure  4  gives  the  additional  factor  to  be  used  when  the  upward 
slant  path  spectral  amplitude  is  desired.  It  is  evident  that  the  low 
frequencies  are  accentuated  and  the  intermediate  frequencies  are 
attenuated  for  an  upward  slant  path. 

Figure  5  shows  the  h  factor  for  a  downward  slant  path.  Again 
the  low  frequencies  are  accentuated  and  the  intermediate  and  low 
frequencies  are  attenuated. 


15 


DOWNWARD  Path 


Fig.  5.  The  additional  factor  for  downward  slant 
path  phase  difference  spectra. 
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Ihe  accentuation  ot  tne  low  Trequercies  is  aue  to  uit  

of  outer  scale  effects  at  the  higher  altitudes,  and  the  attenuation 
of  the  intermediate  and  high  frequencies  is  due  to  the  smaller  turbu¬ 
lence  strength  at  the  higher  altitudes. 

The  horizontal  path  differential  path  contribution  is  essentially 
constant  (3432-2)  along  the  path  for  high  temporal  frequencies  and 
hence  the  upward  and  downward  path  high  frequency,  slant  path  factors 
are  the  same  as  can  be  seen  by  comparing  Figs.  1  and  2. 

To  summarize,  temporal  phase  difference  calculations  have  been 
extended  to  include  slant  and  vertical  upward  and  downward  paths. 

The  calculations  show  that  spectra  are  greatly  influenced  by  the 
variation  of  outer  scale  and  turbulence  strength  along  the  path. 


SUMMARY  AND  CONCLUSIONS 

This  report  covers  work  performed  on  Contract  F30602-72-C-0305 
between  July  and  October,  1972.  The  work  performed  is  in  support  of 
the  atmospheric  propagation  program  at  RADC.  Two  items  are  discussed. 
The  first,  a  theoretical  description  of  the  device  used  to  measure 
angle  of  arrival  variance  at  RADC  is  the  conclusion  of  a  project 
started  son.i  time  ago.  The  analysis  was  performed  to  obtain  a  better 
understanding  of  the  apparatus  and  its  performance.  Results  show  that 
the  device  measures  a  quantity  different  from,  but  with  symmetry 
properties  similar  to,  those  commonly  defined  as  angle  of  arrival. 
Further,  only  light  entering  the  aperture  within  a  small  distance  of 
che  rim  contributes  to  the  fugle  of  arrival  measurement. 

The  second  item  is  temporal  spectra  of  phase  difference 
fluctuations.  This  is  an  extension  of  previous  work  and  covers  slant 
and  vertical  paths.  The  results  show  that  the  temporal  spectra  are 
strongly  influenced  by  the  variation  of  C^  and  L0  along  the  path,  and 
that  the  effects  are  different  for  upward  and  downward  paths. 
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